Collective Enhancement and Suppression of Excitation Decay in Optical Lattices 
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We calculate radiative lifetimes of collective electronic excitations of atoms in an infinite one di- 
mensional lattice. The translational symmetry along the lattice restricts the photon wave vector 
component parallel to the lattice to the exciton wave number and thus the possible emission direc- 
tions. The resulting radiation damping rate and emission pattern of the exciton strongly deviates 
from independent atom. For some wave numbers and polarizations the excitons superradiantly decay 
very fast, while other excitons show zero radiation damping rate and form propagating meta-stable 
excitations. Such states could be directly coupled via tailored evanescent fields from a nearby fiber. 
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Ultracold atoms in optical lattices are now one of the 
most nourishing directions of experimental and theoret- 
ical quantum physics [ij. Beside fundamental physics 
issues on entanglement, measurement and decoherence, 
more and more applications and connections appear as 
test-systems for a wide range of puzzling effects in solid 
state physics. The trapped atoms can be considered as 
artificial crystals, which differ from solid crystals in their 
precise and easy controllability of occupation number, 
lattice constant, lattice depth and symmetry. The lattice 
properties can be controlled through the laser intensity, 
wave length and polarizations, and through laser cooling 
the motion can be cooled and confined to the lowest state 
in each potential well corresponding to the first Bloch 
band. The atomic dynamics is generally well described by 
a Bose- Hubbard model @ , which includes atom hopping 
among nearest neighbor sites and the on-site repulsive 
atom-atom interactions. As central prediction one gets a 
quantum phase transition from a superfluid into a Mott 
insulator phase 0, Q with a fixed number of atoms per 
site for a deep lattice. Other condensed matter phenom- 
ena appearing in more complex situations can be studied 
in generalized optical lattice setups as well Q. 

Optical properties of crystal solids are a well studied 
subject. One of the phenomena that strongly charac- 
terize the optical properties of solids is the formation of 
collective electronic excitations (excitons). In molecular 
crystals a local electronic excitation delocalizes among 
the lattice molecules due to electrostatic interactions. 
The corresponding delocalized eigenstates are Frenkel ex- 
citons [f| , which are quasi-particles that propagate in the 
lattice. In previous work we showed that Frenkel likc- 
excitons also exist for ultracold atom lattices in the Mott 
insulator phase [H, 0] • Enclosing the light field within a 
cavity, such excitons and photons with the same wave 
vectors coherently couple to form cavity polaritons @, Q . 

In the present work we study the radiative properties of 
such excitons, which are formed due to resonant dipole- 
dipole interactions. Each exciton is represented by a wave 
that propagates in the lattice with a given wave number 
6] and coupled to the free space radiation field modes 
into which it can decay. For the moment we restrict 
ourselves to an infinite one dimensional lattice with one 



atom per site. This allows one to directly calculate the far 
field radiation pattern and the resulting effective damp- 
ing rate of the exciton explicitly. As a result of the lattice 
symmetry an exciton with a fixed wave number can emit 
a free photon with the same wave number component 
along the lattice axis. In particular we investigate, how 
the damping rate, which significantly differs from the in- 
dependent atom result, changes with the wave number 
and polarization direction. 

We consider a one-dimensional string of two-level 
atoms with an electronic transition energy E a = fku a 
with lattice constant a and one atom per site (see fig- 
ure (1)). Energy is transferred among atoms at different 
sites by resonant dipole-dipole interactions @|, so that 
the electronic excitation Hamiltonian contains two terms, 
the on-site excitation and the energy transfer, it reads 

H ex = BlB n + Mnm{0) B^ n B m , (1) 

n nm 

where B^ n and B n are the electronic excitation creation 
and annihilation operators at site n, respectively, which 
are assumed to be boson operators at low excitation den- 
sity with the commutation relation [B n , B^] = S nm . The 
Hamiltonian can be easily diagonalized, in exploiting the 
lattice symmetry, by using the transformation 

B « = ^I> 4fc ^' ( 2 ) 

where N is the number of sites, and z n = an is the posi- 
tion of site n. We get the collective electronic excitation 
Hamiltonian 

H ex = J2^ex(k,e) B\B k , (3) 

k 

with the dispersion uj ex (k,9) — u a + J2l J{L, 9)e lhL , 
and L = (n — m)a. Restricting the energy trans- 
fer to only nearest neighbor sites, we explicitly obtain 
w ex (k,9) — u) a + 2 J (9) cos(ka), where we used J(9) = 
J (a, 9), and for the dipole-dipole interaction we have 

HJ{9) — 4 J^ oa i (l — 3 cos 2 6>), here fi is the electric tran- 
sition dipole, and 9 is the angle between the transition 
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dipole and the lattice axis, as seen in figure (1). Includ- 
ing corrections from more distant neighboring atoms will 
change this spectrum and the particular form of the ex- 
citons, which cannot be calculated explicitly any more. 
Nevertheless, the results and conclusions of the following 
life time calculations are generally only weakly influenced 
by these corrections as the wave nature of the exciton is 
mainly determined by symmetry properties of the lattice. 
In using periodic boundary condition, the wave number k 
takes the values k = jfcp, with p = 0, ±1, ±2, • • • , ±N/2. 
It is clear that in place of the discrete atomic transition 
we get an energy band of band width 4J(9). 

Due to the energy transfer, an quasi-stationary elec- 
tronic excitation is dclocalizcd in the lattice and best rep- 
resented by a wave that propagates to the left or the right 
direction with wave number k, which is a good quantum 
number. Such collective electronic excitations are called 
"excitons". While some excitons directly couple to in- 
coming plane waves, for other excitons with a fixed wave 
number k and a prescribed transition dipole direction 9, 
such a straightforward excitation is not possible. As has 
been implemented recently, one way to overcome this is- 
sue is the use of a thin pulled fiber along the lattice @. 
The atomic lattice is generated parallel to the fiber at 
given distance, where the fiber photons and the corre- 
sponding lattice excitons are coupled. Hence, by sending 
photons with prescribed wave number and polarization 
through the fiber we can excite an exciton in the lattice 
with the same wave number and polarization through the 
evanescent wave coupling [Io| . 
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FIG. 1: A one-dimensional lattice of lattice constant a, and 
with one atom per site. The transition dipole fi is plotted 
which makes an angle 6 with the lattice axis. 

Here we calculate the electric field and intensity of the 
radiation emitted by a ID lattice exciton into free space. 
The free-space radiation field Hamiltonian is given by 



Hr 



qA 



(4) 



where b ^ and bq\ are the creation and annihilation op- 
erators of a photon with wave vector q and polarization 
A, respectively. The photon energy is E p h(q) = huj r (q) = 
hcq. The electric field operator is 



E ra d(r) = i 



qA 



bq\ e qA e zq r - b 



t * 

'qA e qA« 



**}, 

where e q > is the photon polarization unit vector, and V 
is the normalization volume. 

The light-matter coupling is given by the electric dipole 
interaction Hi = —p, ■ E ro( j, where the atomic transition 



dipole operator is defined by ft — ftJ2n {^n + B^). We 
separate the photon wave vector into two components, 
one parallel to the lattice axis and the other orthogonal 
to it, namely we use q = q' + q z z, with q 2 = q' 2 + q 2 z . 
The atomic site positions are given by r„ = z„z. We 
represent the electronic excitations by excitons, in using 
the inverse transformation of equation @ , and by using 
the lattice symmetry property -k J2 n e*^~ fc ) z ™ = 5 Qzk . 
The coupling between the excitons and the free radiation 
field in the rotating wave approximation and with linear 
polarizations, is given by 



H, 



^2 ih 9^ { b i' 



}• 



where now q = (q', k), with the coupling parameter 



% q A = - 



hcqN 
2^V 



iP> ' e qA) ■ 



(6) 



(7) 



It is seen that the coupling is only between excitons and 
photons with the same z component wave vectors. The 
translational symmetry along the lattice axis results in 
a conservation of momentum along the axis. Hence an 
exciton with wave number k couples only to photons with 
the same k parallel to the lattice. Such a result strongly 
affect both the radiation emitted by an exciton with a 
fixed k, and the damping rate of such exciton, as we show 
in the following. The effect of similar retarded interaction 
on the exciton spectrum at low dimensional molecular 
crystals was studied by Agranovich et. al. [ill ]. 

Now we calculate the electric field and the intensity 
of the light emitted from one dimensional optical lattice 
of cold atoms. We start from the radiation field opera- 
tor equation of motion, solving it formally and substitut- 
ing the source part back in the electric field operator of 
equation ([5]), then we obtain for the positive part of the 
electric field, for a fixed k and 9, 



q'A 

t 



u} r (q)VN ei [q'.p + kz^Ak)t] 
2e Q V 



e q A (M -e q A) / dt' B k (t') e *KM-«v to] (*-*') , (8) 



where the exciton operator in a rotating frame is defined 
by B k (t) — B k (t) e - lw e( fe )* 5 an d the observation point is 
at r = p + zz. As the wave vector component parallel 
to the lattice is fixed, we have only to sum over wave 
vectors in the plane normal to the lattice direction. The 
sum over q' casts into the integral 



E 



s 

4^2 



d 2 q' = 



S 
4^2 



d<\> \ q'dq', (9) 



where V = SL, and S is the normalization plane, with 
the lattice length L — Na. We use also the summa- 
tion over the photon polarization X^A e qA e qA = 1 — 
For the different vectors we choose: the wave vector is 
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q = (q' cos 4>, q' sin <fr, k), the transition dipole is ft = 
(/i sin#, 0, /i cos#), and the observation point is chosen 
to be at r = (p, 0, z), as seen in figure (2). 

After doing the angular integral we keep only the far 
field terms proportional to 1/ ^fp. As solvable example 
we consider long wave length (small wave number) exci- 
tons with Then we apply the Weisskopf-Wigner 

approximation, which is in the spirit of the Markov ap- 
proximation [12j. The retarded electric field is given by 



E+(r, *) = (!+<; 



ineoac 3 / 2 \ Np 
and the intensity operator is 



(10) 



pcos9\ 2u) e (k) 3 
4e a J Nirpc 3 " k 



B{(t-p/c)B k (t-p/c). 



(11) 



The expectation values of the exciton operator is 
(B k {t-p/c)) = (.B fe (G)) e -^e(fc)(*-p/c) e -r fc (t- P /c)/2 j 

where r k is the damping rate of an exciton with wave 
number k. The expectation value of the intensity is 



<*> = 



^cos#\ 2oj e (k) 



4e n a 



Nnpc 



.. ,4(0)^(0)) e' T ^-Pl c \ 



(13) 

Our next task should be to calculate the exciton damping 
rate, Tk, into free space, and to compare the result with 
the single atom damping rate in order to emphasize the 
cooperative effect among the atoms on the damping rate. 



X 


(p,0,z) 

i q 


V 

\a j 













FIG. 2: The lattice is in the z direction, where the transition 
dipole p is in the (x — z) plane and makes and angle 6 with 
the lattice direction. The exciton has a wave number k along 
the lattice. The photon wave vector q has a component q' 
normal to the lattice in the (a; — y) plane, and a z component 
k parallel to the lattice which is equal to the exciton wave 
number. The observation point is at (p,0,z). 

The radiative exciton damping rate with a fixed k can 
be evaluated using the Fermi Golden Rule, which reads 



2tt 



Y,\(f\Hi\i)\ 2 5(E ex (k)-E ph (q)) 

q'A 



(14) 



The initial state is of a single exciton of wave number 
k in the lattice and an empty radiation field, that is 



\i) = \lex(k),0 p h), and the final state is of zero exci- 
tons and a single radiation field photon of wave vector 
q, that is |/) = \0 ex , l p ^(q)). The matrix element is 

given by = (p ■ e qA ). As before the 

emitted photon has a wave vector component parallel 
to the lattice equals to the exciton wave number k. The 
summation is over the wave vectors q' normal to the 
lattice direction. The delta function cares for the con- 
servation of energy of the exciton and the emitted pho- 
ton. The summation over the photon polarizations is 

obtained by the relation ' e qA| 2 = Im| 2 — q ^ ■ 

The sum over q' is converted into the integral of equa- 
tion ([9]). As before we assume a transition dipole with 
a fixed direction of p, = (p sin#, 0, p, cos 9), and with 
q • pZ = q'psm9cos(f> + kp cos 9. After integration we 
obtain the result 



k?EUk) 



1 



cos 



(tick) 2 



(2 cos 2 9 - sin 2 



(15) 

Note that the free atom damping rate is r at = 3 ^ egti d c 3 ■ 
As expected for long enough wave vectors k and cer- 
tain angles 9 the damping rate can be much larger than 
the single atom one, as we get a superradiant enhance- 
ment analogous to the Dicke model. However, we see 
that for some k and 9 the damping rate can be much 
smaller than the single atom one or even completely van- 
ish. Such metastable states (dark states), which appear 
for lattices with two atoms per site [8], thus also exist 
in infinite chains of atoms as treated. Quite generally at 
9 — 0° we can calculate a critical wave vector k c , given 
by E ex (k c ) = hck c , above which the damping rate is zero 
and the excitations can no longer decay radiatively. For 
general polarization such a critical wave vector exists if 
the following equation has a solution 



(hck c 



1 + cos 2 9 



E 2 x (k c ) 2 cos 2 9 - sin 2 9 ' 



(16) 




FIG. 3: The damping rate vs. ka, for (a) 9 — 0, and (b) 
6 = 90. The dashed-line is for a single atom damping rate. 

Let us now exhibit this behavior in several plots of the 
damping rate. We use the parameters: the transition 
energy is E a = 1 eV, the lattice constant is a = 1000 A, 
and the transition dipole is fJ, = 1 eA. In figure (3. a) 
we plot the damping rate as a function of ka for 9 = 0°, 
the plot includes also the damping rate of a single atom 
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(dashed line). Note that for small wave vectors, ka ~ 
the exciton damping rate is much larger than the free 
atom damping rate, and these states are superradiant 
states. With increasing wave vectors the damping rate 
drops and reaches the single atom rate, and beyond a 
critical wave vector k c the damping rate becomes zero 
and no damping is obtained beyond k c . In figure (3.b) 
the plot is for 9 = 90°, here the damping rates are larger 
than the single atom one for small wave vectors, ka ~ 
0, but now the damping rate increases with increasing 
the wave vector. Namely, the large wave vector states 
became more superradiant states. Exist an angle, around 
9 = 54.7°, where the damping rate changes its behavior. 
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FIG. 4: The damping rate vs. 6 (full-line), for ka = 0.01. 
The dashed-line is for a single atom damping rate. 
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FIG. 5: The damping rate vs. 9, for (a) ka — 0.5, and (b) 
ka — 1. The dashed-line is for a single atom damping rate. 

Next for different values of ka we plot the damping 
rate as a function of 9. In figure (4) we use ka = 0.01, 
here the damping rate is larger than a single atom one 
for all angles, but around 9 = 90° the damping rate is 



lower than that of 9 — 0°. In figure (5. a) we use ka = 0.5, 
now a dramatic change is obtained, where the damping 
rate is smaller than the single atom one around 9 ~ 0° 
and increases to a maximum at 9 — 90°, that is much 
larger than the single atom case and which is strongly 
superradiant state. For larger ka, in figure (5.b) we use 
ka = 1, now for angles up to 45° the damping rate is 
zero, and they are nonzero for 9 ~ 45° — 135° with a high 
maximum at 9 — 90°. For ka — tt at the boundary of 
the Brillouin zone, the non-zero region is for 9 ~ 54.7° — 
125.3°, and with higher rate at 9 = 90° than the previous 
case of figure (5.b). 

We showed that the damping rate and the radia- 
tion emitted of an exciton in one dimensional lattice is 
strongly different from independent atoms. The excitons 
formed by the energy transfer among the lattice atoms 
due to dipole-dipole interactions can be characterized by 
their wave numbers along the lattice. This restricts the 
spontaneously emitted photons to the same wave vector 
component parallel to the lattice. The radiative damp- 
ing rate of such excitons into free space strongly depends 
on both, the wave number and the polarization direc- 
tion of the exciton. In close analogy to superfluorescent 
scattering, some wave numbers and polarizations exhibit 
a much larger damping rate than an independent atom 
(superradiant excitons), while for other wave numbers 
and polarizations the damping rate is much smaller and 
even can vanish completely, so that the excitons now are 
metastable. The present results can be adopted to any 
periodic distribution of optically active materials, e.g. a 
lattice of quantum dots, one dimensional molecular crys- 
tals, periodic semiconductor nanostructures, or a chain of 
trapped ions, etc. Interestingly such excitons might still 
be able to couple to and decay into a nearby material 
structure if it supports the corresponding wave vectors. 
As a single atom can only bear one excitation, several of 
these excitons will interact and might open possibilities 
of nonlinear light interactions at very low intensities. 
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